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Abstract 

Shapere and Wilczek recently found some singular Lagrangian systems which 
O-j spontaneously breaks time translation symmetry. The common feature of their 

models is that the energy functions are multivalued in terms of the canonical 
phase space variables and the symmetry breaking ground states are all located 
at the brunching point singularities. By enlarging the phase space and making 
use of Dirac's theory on constrained Hamiltonian systems, we present the Hamil- 
tonian description of some of the models discussed by Shapere and Wilczek and 
found that both the multivaluedness and the brunching point singularities can 
be avoided, while the spontaneous breaking of time translation becomes more 
transparent. It is also shown that the breaking of time translation is always 
accompanied by the breaking of time reversal. 

Keywords: spontaneous breaking of symmetry, time translation, time reversal, 
constrained Hamiltonian systems 



PACS: 45.20.Jj, 11.30.Qc 

1 Introduction 



Conventional wisdom says that every classical conservative mechanical system has a 
conserved total energy H which generates a continuous time translation symmetry. 
Until very recently, it is commonly believed that time translation symmetry cannot 
be spontaneously broken. This picture has changed drastically due to the work [1] 
of Shapere and Wilczek, who found that in some special Lagrangian systems the time 
translation symmetry can be spontaneously broken in the lowest energy state (or called 
ground state). All models considered by Shapere and Wilczek bears the same property 
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that the energy functions cannot be made singlevalued in terms of the generalized co- 
ordinates and their conjugate momenta. Historically, such systems are called singular 
Lagrangian systems [2, 3] or multivalued Hamiltonian systems [4], and, even though a 
plausible quantization procedure exists [4] in certain special systematic Hamil- 

tonian description for the classical case is still missing. 

In this paper, we shall address the Hamiltonian description of singular Lagrangian 
systems with spontaneously broken time translation symmetry. Before doing this, 
we would like to pay some words clarifying the meaning of spontaneous breaking of 
time translation symmetry. First of all, in Hamiltonian mechanics, states of motion 
are represented by points in the phase space. With properly chosen coordinates on 
the phase space, we can say that a state of motion is just a point in phase space with 
definite coordinate values. The physical notion of motion is just the change of the state 
in the course of time governed by the classical equations of motion, and the ground state 
is the state that minimizes the Hamiltonian function, usually subjects to no motion. 
For standard canonical Hamiltonian systems, the most commonly used coordinates 
on the phase space are generalized coordinates and their conjugate momenta, which 
constitute a set of Darboux coordinates on the phase space regarded as a symplectic 
manifold. However, this does not mean that Darboux coordinates are the only viable 
coordinates on the phase space. In principle, every set of 2n (n being the number 
of degrees of freedom of the mechanical system) linearly independent functions of the 
Darboux coordinates can be chosen as a coordinate system on the phase space. Let 
{Xi,i = 1 • ■ ■ 2n} be any coordinate system on the phase space, under which the 
ground state has coordinates {Xf }. A state is subject to nontrivial motion if there 
exist some i such that = -^Xi ^ 0. Recall that under infinitesimal time translation 
t t + e, Xt(t) -> Xi(t + e) ~ Xi(t) + eXi(t). If Xf = for all i, then the ground 
state is preserved by time translation symmetry. This is the case for most familiar 
mechanical systems (e.g. classical harmonic oscillator). If, on the contrary, Xf ^ 
for some i, then the ground state is not preserved by time translation, i.e. the time 
translation symmetry is spontaneously broken. We see that spontaneous breaking of 
time translation symmetry requires nontrivial motion in the ground state. 

One may wonder how a true ground state of classical conservative mechanical sys- 
tem can carry through nontrivial motion. This question is beautifully answered by 
Shapere and Wilczek by presenting concrete examples. The models they discussed are 
all Lagrangian systems depending on higher powers of generalized velocities. For such 
systems the conjugate momenta are nonlinear functions in terms of generalized veloci- 
ties, and this makes the Hamiltonian multivalued in the conjugate momenta. The trick 
is that the ground states always appear at the branching points (they called these the 
"turning points") so that their motion is simultaneously governed by several different 
Hamiltonian functions, rendering the result of the motion unpredictable. From the 
point of view of Hamiltonian mechanics, the multivaluedness of the Hamiltonian sig- 
nifies that the system is ill defined, and the ground states at the branching points of 
the multivalued Hamiltonian are in contradiction with the usual smooth minimization 
procedure. Below we shall present Hamiltonian descriptions of some of the models con- 
sidered by Shapere and Wilczek, avoiding the appearance of multivalued Hamiltonian 
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and the branching point singularities. The only price needed to be paid is that we have 
to abandon the use of Darboux coordinates on the phase space. 



2 General considerations 



Consider the following Lagrangian involving higher powers of generalized velocity, 

n 

fc=l 

where denotes the generalized coordinate on the configuration space and fk{4>) are 
functions of with / n (0) > and at least one of the fk{4>) is negative for k = 1, n— 1. 
This Lagrangian contains only even powers of because we implicitly imposed a time 
reversal symmetry. 

The conjugate momenta associated with reads 

5 <t> k =i 

Following the standard procedure, we introduce the Hamiltonian 

*L ok — 1 

n = w> - c = £ -^j-f^ 2 * + fM- 
k=i 

The Hamiltonian should in principle be regarded as a function of the phase space 
variables ((f), p)- However, since p is a nonlinear polynomial in <fi, 7-L will in general be 
multivalued in p. So, for a given set of values of (4>,p), there can be several different 
values for H. If, as the states tend to some special points (<p^ B \p^), some of the 
different values of % become degenerate, we call these special points the branching 
points of the Hamiltonian. 

The existence of branching points indicates that the original Lagrangian system is 
singular and possesses some special features. Indeed, though at the first glance, the 
Lagrangian depends only on (0, 0) and the the standard Euler-Lagrangian equation of 
motion 

5C d (bC 



5<p dt V^0_ 

should hold without problem, a moment thinking would reveal that things are more 
involved. Consider the action S = J Cdt for the above system. After integration by 
parts, the Lagrangian can be changed into a form containing higher order time deriva- 
tives of 0. In the conventional treatment of such systems, one takes 0, 0, ^ e _ 1 
as independent generalized coordinates and define variational derivatives of the cor- 
responding Lagrangian with respect to the time derivatives of these coordinates, i.e. 
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0, 0, 4jj(f), as generalized momenta {£ is the highest order of time derivatives appear- 
ing in the Lagrangian) . In doing so, one effectively turns the system into a constrained 
Lagrangian system with a larger number of degrees of freedom, and to get the consistent 
time evolution of the unconstrained degrees of freedom, Dirac's theory for constrained 
Hamiltonian systems [5] is often employed. 

Actually, the idea of turning the Lagrangian system into a constrained Hamiltonian 
system is the key to resolve the problem of presenting a consistent Hamiltonian de- 
scription for singular Lagrangian systems with spontaneously broken time translation 
symmetry. To employ Dirac's method, it is not necessary to carry out the integration 
by parts in the action. Rather, we can manually introduce novel coordinates 7, p and 
introduce a new Lagrangian 

n 1 

1 = E Yk m)p2k - /o(0) + 7(p " ^ (2) 

fc=l 

which is equivalent to the original one by solving the 7 equation of motion and sub- 
stituting back in L. Since the Lagrangian (2) depends only linearly on and does 
not depend on 7, p, it is clear that this system is a constrained system with a larger 
phase space containing the phase space of the original system as a subspace. Our task 
is then turned into finding the consistent Hamiltonian structure on an unconstrained 
subspace of the enlarged phase space which gives the proper Hamiltonian description 
of our problem at hand. This last step, however, must be carried out on a case by case 
basis. 

Let us remark that although in the above we considered only systems involving 
a single generalized coordinate, the generalization to cases with multiple generalized 
coordinates is straightforward. 

3 Model analysis 

To make our analysis clear and understandable, we shall consider some of the models 
presented by Shapere and Wilczek, i.e. the 4 model, the fgh model and the double 
sombrero model. These models will be analyzed in separate subsections. 

3.1 The 4 model 

The 4 model is unambiguously defined by its Lagrangian 

m = + 

where A, k > are constant parameters. Clearly, this is the simplest case of n = 2 in 
(1) with the choice f'2 = X, fi = —k and f = 0. Using the idea outlined in the last 
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section, we introduce the following Lagrangian involving two novel coordinates p and 

7, 



L{(t>, P, 7) = — + 4 p4 + 7 ^ p _ ^ 

This system possesses an enlarged phase space spanned by (p, 7, (ft, ir p , 7r 7 , 7r^), where 
(7T p , 7r 7 , 7r^) are conjugate momenta of the respective generalized coordinates. This en- 
larged, 6-dimensional phase space does not possess a well defined symplectic structure, 
because the definition of the conjugate momenta gives rise to 3 primary constraints, 

Gi = Titf, + 7 ~ 0, 
G 2 = vr 7 ~ 0, 
G 3 = ti p ~ 0. 

Using standard Legendre transformation we can obtain the expression for the Hamil- 
tonian H associated with L(<p, p, 7). Consequently, we are ready to analyze the time 
evolution of the primary constraints with respect to the total Hamiltonian 

3 

H total = H + ^2 ViGi 

i=l 

K A 

= -p 2 - ^p 4 - IP + faGi + p 2 G 2 + ^3^3 (3) 

using the canonical Poisson brackets for the phase space variables (p, 7, 0, n p , 7r 7 , 7r^). 
It turns out that only the time evolution of G3 leads to a secondary constraints 

G 4 = -up + Xp 3 + 7 ~ 0, 

and there is no further constraints. Notice that the Lagrangian multipliers pi and jli 
appearing on the first and second lines of (3) are not the same quantities. We have 
absorbed the unknown 0, 7, p etc into the variables pi. 

Next we need to evaluate the matrix Poisson brackets between the constraints, i.e. 

{G a , Gp}. 

\ 

-1 
k- 3Xp 2 ' 

k + 3Ap 2 / 

As long as k — 3Ap 2 7^ 0, the above matrix has rank 4 and all the constraints G\ ~ G4 
are of the second class. Following the standard procedure we introduce the Dirac 
Poisson bracket 

{A, B} DB = {A, B} - {A, G a }(M~ 1 ) al3 {G j 3, B} 



M a(3 = 

The result reads 



(0 


1 


-1 











^0 





1 - 
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for any two functions A, B on the phase space. After doing so, all second class con- 
straints can be consistently set equal to strong zeros. The only remaining phase space 
variables are (0, p), with 

{>, P}dB = 777 q , {P, P}DB = {4>, 4>}dB = 0. (4) 

6Xp z — k 

The two dimensional subspace spanned by (</>, p) will be henceforth referred to as 
reduced phase space. The final Hamiltonian on the reduced phase space is given by 

1 3 

H fin = ~2 K P 2 + 4 V' 

and the Hamiltonian equations of motion are 

<[> = {0, H fin } DB = p, p = {p, Hf in } DB = 0. 

The final Hamiltonian is smooth and single valued, its minima obey the usual extremal 
conditions —q^- = = 0. The minima occur at 

p G = ±v^73A, 

and G is not fixed by the minimization condition. Rather, its rate of change in time 
is given by the equation of motion, 

G = p G = ±vV 3A - 

If we take (0 , p ) as parameters specifying the ground states, then the above system 
will possess two distinct families of infinitely degenerated ground states parametrized 
by (0 G ,p G ) = (v^73At + O , y^/SX) and (0 G ,p G ) = (-^/3\t + O , ~^N>) re- 
spectively, where the constant 0o is the initial value of <p G . The explicit dependence of 
cf) G on t indicates that the time translation symmetry is spontaneously broken. More- 
over, the actual ground state must be chosen from only one of the two families of 
degenerate ground states. The concrete choice of ground state also breaks the time 
reversal symmetry. 

What happens when k — 3Ap 2 = 0? It appears at first glance that in this case 
the matrix M becomes uninvertible and the above procedure for constructing Poisson 
brackets on the reduced phase space breaks down. This is also signified in the final 
Poisson bracket {<f), p} = 3Ap 2_ K - If the singularity at the right hand side appears at 
other places in the phase space, this would indeed be a problem. However, notice that 
at exactly this particular singularity, the final Hamiltonian acquires its critical value. 
In the standard treatment of Hamiltonian mechanics, the symplectic structure on the 
phase space is often written as 

Q = dq A dp = dt A dH, 

where in the second equality the time parameter t and the Hamiltonian H are taken 
as a pair of action-angle variables. It is clear that the the critical points of H, the 
symplectic structure vanishes and the Poisson bracket is not well defined. Nonetheless, 
such singularities will not affect the effectiveness of the Poisson bracket in general, as 
we have witnessed above. In particular, when evaluating the equations of motion using 
the Poisson brackets (4), all such singularities disappears and the resulting equations 
hold without problem even at the criticalities of the Hamiltonian. 
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3.2 The fgh model 



The previous model does not allow the inclusion of a potential, which makes it not 
useful when interaction is involved. To allow the inclusion of interacting potentials, let 
us look at a slightly more complex model - the fgh model specified by the Lagrangian 

with > 0,g(<p) < 0. This model is the special case of (1) with n = 2, f 2 = 

4/(0), fi = 2g(<f)) and /o = h(<f)). Minimization of the corresponding energy function 
(here and below we shall omit the arguments of /, g and h and denote the derivatives 
such as -h as for simplicity) 

E = 3f<p 4 + g<p 2 -h, 

requires both <fi 2 + ^ = and + h = c, where c is a constant. This is impossible 
for generic h, because the former condition implies nontrivial time dependence of 0, 
while the latter requires to take some specific constant value(s). To circumvent this 
contradiction, it is necessary [1] to take h to be the solution of the algebraic equation 

2 

jTjj + h = c. Doing so the Lagrangian is rewritten as 

L{ ( j ) ) = f^ + g^ + c-^-. 



We proceed as outlined in the last subsection and introduce an equivalent con- 
strained Lagrangian, 



L{<f>, P, l) = fp A + 9P 2 + c- — + j{p - 0). 

After a standard procedure for constructing Dirac Poisson bracket, we are left with a 
reduced phase space spanned by (p, <fi) with the following Poisson bracket and Hamil- 
tonian, 

{</>, p}db = 2 ^ 2 - ^ , {0, 4>} = {p, p} = 0, (5) 

4 , 2 , 9 _ f ( 2 , 9 



H fm = 3/p 4 + gp 2 - c + ^ = 3/ \^p 2 + f- f ) (6) 



The Hamiltonian equations of motion are 

(j) = {<j),Hf in } DB = p, 



I i 

P = {P> H fin } DB = (f^g - 2g^f) - —f 
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Once again, the apparent singularities in the {0, p}db bracket appears at the critical 
points go Hfi n , which is not a problem as we have explained in the last subsection. 

Eq.(6) shows that the final Hamiltonian is smooth and single valued, and its minima 
are not isolated points, but rather a pair of curves in the reduced phase space given by 

the equations p = ±y^. It is easy to check that each of these two curves is consistent 

with the Hamiltonian equation of motion. 

Due to the complicated form of the equations of motion, it is not easy to determine 
the time dependence of <p G , p G in the ground states. However, given any value of 

(p G consistent with the equation of motion, the states (0 G ,p G ) = (^(f) G , zby ^^^ j 

constitute two families of infinitely degenerated ground states. The time dependence 
of each states in these two families indicates that the time translation symmetry is 
spontaneously broken. Moreover, the actual ground state must be chosen from only 
one of the two families. The concrete choice of ground state automatically breaks the 
time reversal symmetry. 



3.3 The double sombrero model 

The double sombrero model is a little bit more complicated than the previous two 
models because more degrees of freedom are involved. The model Lagrangian is 

l( p a) = \(p 2 + p 2 <P 2 -^ + ^p 2 -^p\ (7) 

where k, /i, A are all positive constants. Following the above procedure, we first intro- 
duce the first order Lagrangian, 

L(p, (f), tp, ip, 7, 77) = ^(ip 2 + p 2 ^ 2 - k) 2 + ^p 2 - jp 4 + j((p - p) + rj(ip - 

The construction of Dirac Poisson brackets in this model is much involved comparing to 
the case of the previous two models because the presence of more constraints. However, 
with the help of existing algorithm [6], it is not difficult to write down a computer 
algebra procedure using Maple. It turns out that the reduced phase subspace is 
spanned by (p,<p,(p,ip), with the following Dirac Poisson brackets: 

0, 

3p 2 ^ 2 + <p 2 - k 
(3p 2 ip 2 + 3<£> 2 — n)(p 2 ijj 2 + ip 2 — k) ' 

(3p 2 ip 2 + 3y2 2 — K){p 2 ip 2 + Lp 2 — 

(3p 2 ip 2 + 3(f 2 — K,)(p 2 ip 2 + if 2 — n) ' 

3p 2 jj 2 -j- if 2 — K 

(3p 2 ip 2 + 3<p 2 — n)(p 2 ip 2 + ip 2 — k)p 2 ' 



{P,<P\db = 
{p,¥>}db = 

{p,iP}db = 

{0, <p}db = 

{0, i'}DB = 



t n 2^ 2p 2 i; 2 + <p 2 -K 

{<P,iP}db 



p (3p 2 ijj 2 + 3<p 2 — k)(p 2 i(j 2 + (p 2 — k) ' 
and 

{p,p}db = {<PA}db = W,¥}db = {iP,iP}db = 0. 
The final Hamiltonian on the reduced phase space reads 

H fm = ^ (3pV 2 + 3^-^ + 1- (Xp 2 - pf - £ - (8) 

which yields the Hamiltonian equations of motion 

P = {/°> H fm } DB = <p, (9) 
= {0, H fin } DB = ij;, (10) 

V - H fm \ DB - -— + {3p ^ 2 + 3Lp2 _ K)ip2 ^ + Lp 2_ K y ( 12 ) 

Eq. (8) shows that the final Hamiltonian is smooth and single valued. The minima 
of the Hamiltonian are determined by the conditions 

3pV + 3<p 2 - k = 0, (13) 
\p 2 - p = (14) 

together with the equations of motion. Concretely, from (14), we can get 

which are constants. Therefore, from the equation of motion p = <p, we get 

ip G = 0. 



Inserting these results into (13), we get 

= +A l 

3p 



i Xk 



and further, from the equation <fi = ip, we get 



3p 



where 4>n is a constant. 
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Finally, the ground states are identified with the following curves, 

(0 G ,p G ,^ G ,<^ G ) = f^^^ + ^ei^ea^^.oj , 

where t\ = ±1, e 2 = ±1 and are independent of each other. The number of independent 
choices of (£1,62) is 4. Therefore, there are totally 4 distinct families of infinitely 
degenerate ground states in this model, each family subjects to nontrivial motion, 
which spontaneously breaks time translation. The actual ground states must be chosen 
from one of the 4 families. Any such choice breaks time reversal symmetry. 

Remarks: 

1. As in the previous two models, the apparent singularities on the right hand sides 
of (11) and (12) are located exactly at the critical points of Hfi n , so these do not give 
rise to any problem as explained earlier. 

2. The consistency of Hamiltonian equations of motion with the corresponding 
Euler-Lagrangian equation of motion for the <j) 4 model and fgh model can be easily 
checked. For the double sombrero model, the consistency between Hamiltonian equa- 
tions of motion and the Euler-Lagrangian equations of motion seems more complicated 
to check. We leave this check in the appendix. 

3. The Hamiltonian equations of motion for the double sombrero model, eqs. (11) 
and (12), seems to be singular at the ground states because the indeterminate 0/0 
appears on the right hand side. However, this is completely superficial. Actually, if we 
first collect a common denominator and then substituting in the ground state values 
of the phase space variables, the right hand sides of both (11) and (12) vanishes. 



4 Discussions 



Spontaneous breaking of time translation symmetry is an amazing discovery. The ex- 
istence of such mechanical systems indicates that there is a hidden piece of classical 
mechanics about which our understanding is still very poor. Nonetheless, the impor- 
tance of such models cannot be overestimated. On the one hand, these models provide 
explicit counter examples to the inverse of Noether's first theorem (which roughly says 
that any conserved charge generates a continuous symmetry). On the other hand, sim- 
ilar Largrangian systems appear widely in modern theoretical physics such as the k- 
inflation models in cosmology [7], ghost condensation [8], higher curvature and f(R) 
gravity models etc. 

The original treatment [1] of Shapere and Wilczek on the spontaneous breaking of 
time translation symmetry has depended on the use of standard Darboux coordinates 
on the phase space. Using such variables, the multivaluedness of energy functions and 
the branching point singularities seem to be inevitable. However, our work shows that 
both the multivaluedness of energy functions and the branching point singularities can 
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be avoided, the only sacrifice to be paid is that non Darboux coordinates have to 
be adopted in parametrizing the phase space. The use of non Darboux coordinates 
also resolves the puzzling jumping effects at the ground states described by Shapere 
and Wilczek. Actually, from the detailed model analysis presented in the last section, 
it is clear that there is no such jumping at the ground states using the non Darboux 
coordinates. We believe that the motion with momentum jumping in the ground states 
is just a signature indicating that the momentum, as part of the Darboux coordinates, 
is not a good coordinate for parametrizing the phase space in such models. As a by- 
product, we also revealed that the time reversal symmetry is broken by the choice of 
ground states in such models. 

Is it possible to extract Darboux coordinates from our treatment on the Hamiltonian 
description? The answer is yes, but the original problems will reappear in terms of these 
coordinates. To be more explicit, let us take the simplest 4 model as an example. It 
is easy to see that both (qi,pi) = (up — Xp 3 ,4>) and (^2,^2) = (0, — K p + Ap 3 ) are 
viable choices of Darboux coordinates. Using either choices of Darboux coordinates, 
the Hamiltonian becomes multivalued, either in gi or in p 2 . The multivaluedness of the 
Hamiltonian in terms of Darboux coordinates also indicates that these coordinates are 
not good coordinates on the phase space. On this point, please also notice that p\ (p 2 
reps.) is nonlinear in q x (q 2 reps.), therefore (q±,Pi) (or (q 2 ,P2)) does not parametrize 
the cotangent bundle of the configuration space as apposed to the cases of canonical 
Hamiltonian systems. 

For usual Hamiltonian systems, passing from classical to the quantum description 
is straightforward, one only needs to replace the canonical Poisson bracket {q,p}pB = 
1 by the canonical commutator [q,p] = ih. However, for the models discussed in 
the present work, the Poisson brackets between the "good coordinates" on the phase 
space are all nonlinear, hence the canonical quantization process does not work. The 
systematic quantization of such models remain an open problem 1 . Perhaps the path 
integral quantization for the 4 model [4] can provide some clue for the other models. 

Another related issue is whether the spontaneous breaking of time translation sym- 
metry can be dynamical, i.e. can we construct models which initially preserve time 
translation but eventually evolve into a phase with spontaneous breaking of time trans- 
lation symmetry? A moment thinking seems to indicate that this is possible. Take 
the fgh model for example. If we take f((p) to be a positive constant, g(<f)) to be a 
function such that g((p) > for the initial value of <ft but decreasing as <j) evolves in 
time and eventually becomes negative, then the ground states will transit from the 
initial symmetry preserving phase (<f) G ,p G ) = (0o,O) to the symmetry breaking phase 



(0 G , p G ) = yf) G , ± y 6/(^gj ) • We leave the detailed study of this phenomenon for later 
works. 

1 After the first version of this paper appeared in arXiv, we have noticed that Shapere and Wilczek 
(together with Xiong) have gone further on the issue of quantization [9] [10]. 
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Appendix 



In this appendix we check the consistency between Euler-Lagrangian equations and 
Hamiltonian equations of motion for the double sombrero model. 

Variation of the Lagrangian (7) with respect to yields 

(2pp0 + p 2 0)(p 2 + p 2 2 - k) + 2p 2 0(pp + pp4> 2 + p 2 00) = 0. 

Similarly, variation of (7) with respect to p gives 

p(p 2 + p 2 2 - «) + 2p(pp + pp0 2 + p 2 00) + [(Ap 2 - p)p - p0 2 (p 2 + p 2 2 - «)] = 0. 

We rearrange these two equations in the form 

«i0 + hp + ci = 0, (15) 
a 2 (f) + b 2 p + c 2 = 0, (16) 

where 

at = p 2 (3p 2 <p 2 + p 2 - k), b 1 =2p 2 p<p, ci = 2pp0(2p 2 2 + p 2 - k), 

a 2 = 2p 2 p0, 6 2 = 3p 2 + p 2 2 — k, c 2 = (Ap 2 - p)p - p^ 2 (p 2 2 - p 2 - k), 

and take (15) and (16) as a system of linear algebraic equations for and p. The 
solution is easily found to be 

v aic 2 -a 2 ci 

^ = ^ T' 

a 2 6i - ai0 2 

b x c 2 -b 2 c x 

P = — i r- 18 

a 2 bi - ai6 2 

Inserting the expressions for a±, a 2 , b\, b 2 , ci, c 2 and after some algebraic simplifications, 
we find that (17) and (18) are exactly the same as the results which we can get by 
substituting (9) and (10) into (11) and (12). 
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